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Abstract 

An approximate stochastic model for the topological dynamics of the pe- 
riodic triangular Lorentz gas is constructed. The model, together with an 
extremum principle, is used to find a closed form approximation to the dif- 
fusion coefficient as a function of the lattice spacing. This approximation is 
superior to the popular Machta and Zwanzig result and agrees well with a 
range of numerical estimates. 
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1. Introduction 

The diffusion coefficient is perhaps the simplest example of a transport 
coefficient as it describes the transport of mass in a system. The Lorentz 
gas, originally proposed as a model for the movement of electrons in a crystal 
lattice, comprises a single particle moving in a lattice of fixed scatterers, and 
in the absence of a field the electron diffuses through the lattice. 

A variety of methods have been employed to calculate the diffusion coef- 
ficient of the Lorentz gas. Machta and Zwanzig p] used a Markov hopping 
process to generate an analytical diffusion coefficient approximation. A sim- 
ilar approach has been used more recently for the three-dimensional Lorentz 
gas [2j. As the model used by Machta and Zwanzig is stochastic it will be 
examined in more depth in section [3] Morriss and Rondoni [3] used the 
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periodic orbit expansion method to calculate the diffusion coefficient but 
also performed detailed calculations using the Green-Kubo relations and the 
mean squared displacement. Baranyai, Evans, and Cohen j3] calculated the 
diffusion coefficient based on the Green-Kubo relation although only for a 
limited number of densities. Gaspard and Baras [5] calculated the diffusion 
coefficient using an escape rate formalism which has been extended to other 
systems [6j. Although these results have been known for some time there 
remains much interest in diffusion in billiard systems [7j and the periodic 
Lorentz gas El EH EE] . 

2. Lorentz Gas Parameters 

We consider a triangular lattice of hard disk scatterers with a minimal 
spacing of w between the surfaces of adjacent disks, as shown in Figure [T] 
The wandering particle moves in a straight line in the area outside the scat- 
terers until it has a specular collision with a scatterer. As long as the spacing 
is small enough, w < — 2, the horizon is finite and there are no paths of 
infinite length, so a collision must occur. The path of the wandering particle 
can be described by a symbolic dynamics [12] constructed by assigning a 
symbol to the next (relative) scatterer. If the scatterer is a nearest neigh- 
bour then we can assign that event as a short flight. If the scatterer was a 
next nearest neighbour then that path is a long flight (see Figure [T]). Any 
displacement of the wandering particle in the Lorentz gas with finite horizon 
is composed of a combination of short and long flights. Although this ap- 
proach has the flavour of a periodic orbit expansion, at no stage do we use 
periodic orbits, and the similarity is just in the symbolic dynamics used. In 
the stochastic model we develop the symbolic dynamics to become the state 
space. 

3. The Machta and Zwanzig Result 

Machta and Zwanzig [1] construct their simple analytical estimate for the 
diffusion coefficient by replacing the trajectory of a particle with a random 
walk between triangular regions of the lattice called traps. A diagram of the 
trapping region is given in Figure [2] The probability of a transition from this 
region is calculated by considering the volume of phase space that will leave 
the trap in a time At. This leads to an expression for the mean occupation 
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Figure 1: A section of the periodic triangular Lorentz gas including a typical path for 
the wandering particle beginning from the central scatterer. This path becomes the basis 
for the deterministic model that we introduce later. The two symbol state space used to 
label each segment of trajectory is a relative symbol defined by shifting the trajectory so 
that the segment originates from the central scatterer, and then choosing the symbol by 
the next scatterer it hits. Thus a flight from the central scatterer to a nearest neighbor is 
labelled s for short and flight to a second nearest neighbor is labelled / for long. 
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Figure 2: The Machta-Zwanzig trap indicated by the shaded region between three scat- 
terers in the triangular lattice. In Figure [T] there is a trap between every group of three 
scatterers. 



time of a trap as: 

7T / \/3(w -\ 'IT - ~ 

6w 



(1) 



The diffusion coefficient on a two-dimensional isotropic lattice can be ex- 
pressed as: 

I 2 

where I = (2 + w)/y3 is the distance between the traps. This gives an 
approximation for the diffusion coefficient as: 

DMZ ~ n(V3(w + 2)i-2ny (3) 

This derivation relies on the assumption that the process of transitions 
between traps is Markov. The more collisions that occur in the trap, the 
more information about which hole the particle entered by is lost. The state 
space that Machta and Zwanzig use differs from the state space that will be 
employed here. There are similarities in the approaches but the method that 
is outlined here does not require the same Markov assumption. The Markov 
assumption is only justified at very small w whereas the approximation out- 
lined below should be applicable over a wider range of values. 
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Klages and Dellago [33] and Klages and Korabel [H] have developed 
successive systematic refinements based on explicit correlations of the non- 
Markov events to improve the accuracy of the Machta and Zwanzig result. 



4. A Deterministic Model 

The diffusion coefficient is defined in terms of the linear growth of the 
mean-square displacement with time, and can be found using the Einstein 
relation: 

C =.im^M. (4) 

Here Ar(t) is the displacement of the particle as a function of time. For the 
Lorentz gas the dynamics consists of repetitions of a free-flight at velocity 
v — 1, followed by a collision with a scatterer giving the diffusion coefficient 
as a function of the spacing between the scatterers w. If no infinite length 
flights are possible, the Lorentz gas is said to have a finite horizon and the 
free-flights are of two types; flights between nearest neighour scatterers (short 
flights) and flights between second nearest neighbours (long flights). Any 
physical trajectory can be written as a sequence of short and long flights, 
where each flight has a length r« and a unit vector direction Uj. 

To evaluate the diffusion constant from the Einstein relation we need both 
the length of the trajectory and the time taken to travel along it. The total 
time after N flights can be written as the sum of the times for each flight 8ti 
as 

v^e x^ r i N s (r s ) + NAri) 

t N = V5ti = V- = — 5 

z — ' ^— ' v v 

i=l i=l 

where N s is the number of short flights, Ni is the number of long flights, (r s ) 
is the average length of a short flight and (77) is the average length of a long 
flight. The square displacement Ar 2 (t N ) at time t^ is given by 

N 

Ar 2 (t N ) = (J2ri"i) 2 . (6) 

i=l 

Notice that despite the fact that in the Lorentz gas the flights occur in a 
particular order, the value of both the square displacement and the time do 
not depend on that ordering. Thus we can rearrange the order of the terms in 
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the sum in equation [6]to collect together the short and long flights separately 
as 

/ N s N t 

Ar 2 (t) = I y^TjUi + ^r i u: i 

\i=l i=l 

Ns \ 2 / N t \ 2 / 7V„ \ / Nl 

iV s iV s AT; iV; N a Ni 

i=l t^j j=l i^j i=l i=l 

(7) 

Averaging the square displacement over all initial conditions, assuming 
that rj and u>j are uncorrelated and that the average (a;, • w^) = for all 
z 7^ j, we obtain 



(Ar 2 (t)) = iV s (r 2 ) + iV / (r 2 ) (8) 

The assumption that (u>j -loj) = is only likely to miss correlations between 
subsequent events, especially when both events are short flights. The domi- 
nant term in repeated short flights will be when the two flights are in nearly 
opposite directions, and then (u)i ■ u>j) < 0. Repeated long flights will also 
produce negative correlations but these occurrences are much rarer. 

The probability of a short flight is P s = N s /N and the probability of a 
long flight is Pi = Ni/N so that in the limit N — > oo the diffusion coefficient 
becomes 

m(r s ) + Pi(ri)) { ] 

To evaluate this expression for the diffusion coefficient we need values for 
the averages (r s ), (ri), (r 2 ) and (rf), and the values of the probabilities P s 
and Pi = 1 — P s which will all be functions of the spacing w. The physical 
constraints on the system set lower bounds on r s and r\ (and hence on their 
averages), as they cannot be lower than the minimum separation between 
the relevant scatterers. Thus r s > w and ri > y/3(w + 2) — 2. The values of 
the averages could be found numerically or by explicit integration over the 
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billiard measure, as could the probabilities P s and Pi, but it is less difficult 
to calculate the diffusion coefficient numerically. 

Equation [9] gives an upper bound on the diffusion coefficient as the corre- 
lations that are ignored are mostly negative and in the numerator. To gain 
some quantitative understanding of the accuracy it is perhaps interesting to 
consider the numerical values obtain for various quantities at a single value 
of the spacing, w = 0.2. Here the deterministic model (Equation |9| gives 
D = 0.208 rather than the correct value D = 0.17. This is obtained using 
the average values (r s ) = 0.446, (r,) = 1.88, (r 2 s ) = 0.2576 and (rf ) = 3.56. 



5. Stochastic Model 

Rather than proceeding with the deterministic model of the previous 
section we will use its structure, in particular Equation to construct 
a stochastic model for the deterministic Lorentz gas. The probability space 
for this model will be formed from the phase space of the Lorentz gas together 
with its natural measure. Rather than considering the whole state space we 
will take the Poincare section and develop a discrete model comprised of the 
two states, long flights and short flights. The state space labels transitions 
between iterations on the Poincare surface rather than parameterising the 
position on the surface. The stochastic model constructed in this fashion will 
not be Markovian, as the correlations in the deterministic trajectory persist 
but this is not an issue as only the long run, or stationary, probabilities are 
required. From this it is obvious that all transient behaviour is lost. 

Our model for the diffusion coefficient is given by 

D = (10 ) 
4(P s r s + Pm) 1 UJ 

where the P s and Pi = 1 — P s are the probabilities for short flights and long 
flights. Here r s and 77 are considered to be parameters which may depend 
on the spacing to. There is now no mapping of this stochastic model on to 
random walk on a lattice as the determinism of the previous model has been 
lost. 

As we still have three parameters in the stochastic model, we look for 
a method of reducing this number. To do this we consider looking for an 
extremum of D as a function of the parameters r s and f;. If we minimise D 
with respect to both fi and r s we find that the required value of f; is less 
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than the physical minimum for r h that is y3(w + 2) — 2. So we set f / at this 
lower bound and minimise D with respect to r s alone. 

The calculation of the extremum is straightforward by differentiating 
equation [9] with respect to r s , we find that 

3D = r 2 8 P? + PsPriprs - r t ) 

dr s A(r s P s + riPi) 

which we set to zero. The positive solution of the quadratic in r s gives a 
minimum at 

» h{\fPi-Pi) 



P (12) 

r s 

Notice that r* now depends on w through the dependence of the probabilities 
P s and Pi and the fixed value of rj. The stochastic model now only depends 
upon a single parameter Pi = 1 — P s . After some algebraic manipulation the 
diffusion coefficient is given by 

2P S 2 v ; 

6. Probabilities 

The probabilities of short and long nights should be an integral over the 
regions of initial conditions, weighted with the uniform measure, leading to 
short and long flights and hence are in principal exactly computable. Rather 
than undertake this computation directly, the probability is calculated from 
a numerical simulation which is equivalent to a Monte-Carlo evaluation of the 
integral. We calculate the probability based on the frequency of occurrence 
of short and long flights over a range of spacings w for trajectories of at 
least 5 x 10 5 collisions. The probability of a long flight Pi as a function of w 
appears continuous and smooth and is fitted well by a simple power law in 
w. 

Pi = (3w a (14) 

where the coefficient choices a = 2(a/3 — 1) and = 4g, give a good fit to the 
data. We plot the simulation probabilities and the power law fit in Figure [3] 
Numerical results for the diffusion coefficient are obtained using the fitted 
expression for the probability Pi and the lower bound for f;, thus the final 
result for the diffusion coefficient is 

_ (y/3(w + 2) - 2) (y//?^/ 2 - Pw a ) 

U ~ 2(1-/3^) (15) 
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Figure 3: The probability of a long flight as a function of scatterer spacing w. Each point 
is a molecular dynamics calculation of at least 5 x 10 5 collisions and the solid line is the 
power law fit P; = j3w a . Notice that this fit is a very good approximation throughout the 
whole range of spacings. 
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Figure 4: A comparison of values of the diffusion coefficient obtained from numerical 
simulation with two analytic results, that of Machta and Zwanzig (Equation [3] the dashed 
line), and the result obtained here in Equation f5 (the solid line). The diamonds are Green 
Kubo results from [Tj, the squares are escape rate calculations from Gaspard and Baras 
[5], and the filled circles are Green Kubo calculations from Baranyai, Evans, and Cohen 
[3]. The triangles are the Green Kubo results, and the upside down triangles are periodic 
orbit calculations from Morriss and Rondoni [SJ. 



A comparison of these results with simulation results is given in figure |4j 
This equation fits the other data remarkably well, only falling outside the 
error bars at w — 0.3 and for very small values of w. A deviation may be 
expected at large w as we are using the lower bound for r^. 

At w — 0.2 the stochastic model (Equation 15) gives a value of D = 0.171, 
using r* = 0.343 which less than the value of (r s ) and should therefore lead to 
a smaller value for D. At first glance, the stochastic model appears to replace 
( r s) by ( r s) 2 an d the same replacement for the long flights. The difference 
in these two terms is (r 2 ) — (r s ) 2 = 0.058 and (rf) — (r{) 2 = 0.03 which 
accounts for only part of the difference between the diffusion coefficients of 
the deterministic and stochastic models. 

The fact that the minimum value of the diffusion coefficient gives the best 
value for D over the whole range of w suggests that there is some overall 
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minimum principle at work, although there is no strong justification for this. 



7. Small w limit 

The limiting behaviour of the diffusion coefficient is incorrect in the small 
w limit, where linear behaviour is expected |15| . For small w the power law 



fit for the probability, Equation 14, may be less accurate. The derivation of 
the diffusion coefficient is independent of the form of the probability. For 
the derivation to hold and the diffusion coefficient to be linear as w — > the 
probability must be of the form: 

Pi « aw 2 (16) 

where a is the slope approaching w — 0. Taking the following form for the 
probability we can ind a better fit over the entire range of w. 

P l = /3( e -^ w w 2 + (1 - e^ w )w a ) (17) 

where a and are as given before. The parameter 7 is fitted to the simulation 
data and found as 167.37. With this probability the diffusion coefficient has 
the correct limiting behaviour for small w. 

8. Conclusion 

The simple stochastic model presented here works remarkably well over 
the full range of physical w values. Based on very little information, the 



diffusion coefficient for the Lorentz gas in Equation 15 appears in agreement 
with previous numerical results. The main advantage of the present method is 
computational simplicity and surprising accuracy. A dynamically motivated 
stochastic model, combined with a minimisation procedure and a simple 
power law approximation to the probabilities gives an accurate closed form 
approximation to the diffusion coefficient. It can be corrected at very small 
w to have the correct limiting behaviour. 
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